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The constraint reaction force of ideal nonholonomic constraints in time- dependent mechanics 
on a configuration bundle Q — > R is obtained. Using the vertical extension of Hamiltonian 
formalism to the vertical tangent bundle VQ of Q — > R, the Hamiltonian of a nonholonomic 
constrained system is constructed. 

I. INTRODUCTION 
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This work addresses the geometric theory of nonholonomic constraints in time- dependent 
mechanics. We refer the reader to Refs. 1-7 for the autonomous case. We follow the approach 
based on the D'Alembert principle because the variational methods with Lagrange multipliers 
are not always appropriate to nonholonomic constraints (see Refs. 2,5,6,8). 

Let the jet manifold J l Q be a velocity phase space of time-dependent mechanics on a 
configuration bundle Q —>■ R. The most general nonholonomic constraints considered in 
the literature are given by codistributions S or, accordingly, by distributions Ann (S) on the 
jet manifold J 1 ^. 9 ' 10 Distributions on a configuration space Q and submanifolds of the jet 

^Cj ■ manifold J l Q can be also seen as nonholonomic constraints. 8 ' 11 ' 12 

i > ■ 

Dealing with nonholonomic constraints in time-dependent mechanics, one usually studies 
the following problem. Let £ be a second order dynamic equation on Q and S a codistribution 
on J X Q whose annihilator Ann (S) is treated as a nonholonomic constraint. The goal is to 
find a decomposition, 

Q C = e + r, (1) 

where £ is a second order dynamic equation obeying the condition 

CcAnn(S). (2) 

One can think of £ as describing a mechanical system subject to the nonholonomic constraint 
S, while (— r) is the constraint reaction acceleration. The decomposition (P however is not 
unique. In the case of Newtonian systems, including nondegenerate Lagrangian systems, we 
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obtain the decomposition (|T|) which satisfies the D'Alembert principle for ideal nonholonomic 
constraints. We construct the Hamiltonian counterpart of the constrained equation of motion 
(fj). We show that this can be seen as Hamilton equations in the framework of the vertical 
extension of Hamiltonian formalism to the configuration space VQ which is the vertical 
tangent bundle of Q — > R. This may be a step towards the functional integral formulation 
of nonholonomic time- dependent mechanics and its further quantization. 

II. GEOMETRIC INTERLUDE 

All manifolds throughout the paper are real, finite-dimensional, second-countable (hence, 
paracompact) and connected. 

We refer the reader to Refs. 8-11,13-15 for the geometric formulation of Lagrangian and 
Hamiltonian time-dependent mechanics. In accordance with this formulation, a configura- 
tion space of time-dependent mechanics is an (m + l)-dimensional fiber bundle Q — > R, 
coordinated by {t,q l ). Its base R is treated as a time axis provided with the Cartesian 
coordinate t. With this coordinate, R is equipped with the standard vector field d t and the 
standard 1-form dt. For the sake of convenience, we will also utilize the compact notation 
q x , where q° = t. Obviously, any fiber bundle Q — > R is trivial, but it cannot be canonically 
identified to a product R x M in general. Different trivializations Q = RxM correspond 
to different reference frames. 

The velocity phase space of time-dependent mechanics is the first order jet manifold J X Q 
of Q — > R, coordinated by (t, q\ q\). There is the canonical imbedding, 

A : Jig - TQ, (t, q\ q\) h- (t, q\ i = l,q l = (3) 

of J l Q onto the affine sub-bundle of the tangent bundle TQ of Q which is modelled over the 
vertical tangent bundle VQ of Q — > R. From now on we will identify the jet manifold J l Q 
with its image in TQ. 

Similarly, we have the imbeddings, 

J 2 Q ^ J X J X Q ^ TJ X Q, 

(t, q\ qi q l tt ) ^ (t, q\ qi t = i, q l = ql qt = qi), 

where J 2 Q, coordinated by (g A ,g|,g| t ), is the second order jet manifold of the fiber bundle 
Q — ► R. The affine bundle J 2 Q — ► J X Q is modelled over the vertical tangent bundle, 

VqJ'Q = J l Q x VQ, (4) 

Q 

of the affine jet bundle J X Q — ^ Q. 

The jet manifold J X Q is provided with the canonical tangent-valued form, 

v = Q { <g) dl 
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where 9 l = dq 1 — q\dt are the contact forms. We have the corresponding endomorphism, 



v(dt) = -q\d\, v(di) = 81 v(%) = 0, 

of the tangent bundle TJ l Q and that, 

v(dt) = 0, v(dq l ) = 0, u(dgj) = 6\ 

of the cotangent bundle T*J l Q of J X Q. The nilpotent rule v 2 = holds. 
Due to the imbeddings (^), any connection, 

r = dt®(d t + r%), 

on a fiber bundle Q — > R can be identified with a nowhere vanishing horizontal vector field, 

r = ^ t + r i a i , (5) 

on Q which is the horizontal lift of the standard vector field d t on R by means of V. Con- 
versely, any vector field r on Q such that dt\T = 1 defines a connection on Q — > R. 
Accordingly, the covariant differential, 

D r : JiQ-tVQ, q l oD r = ql-r, 

associated with a connection r on Q — > R, takes its values into the vertical tangent bundle 
of Q -> R. 

Remark: From the physical viewpoint, a connection (j5|) sets a reference frame. There is one- 
to-one correspondence between these connections and the equivalence classes of atlases of 
local constant trivializations of the fiber bundle Q —>■ R, i.e., such that transition functions 
q l — > q n of the corresponding bundle coordinates are independent of t, and T = dt with 
respect to these coordinates. 13-15 In particular, every trivialization of Q defines a complete 
connection r on Q —>■ R, and vice versa. 

A connection £ on the jet bundle J l Q — > R is said to be holonomic if it is a section, 

Z = d t + qid i + edl (6) 
dtj£=l, £J£ = 0, 

of the holonomic sub-bundle J 2 Q — > J X Q of the afline jet bundle J 1 J x (5 — » t/ 1 ^. Holonomic 
connections (^) make up an affine space modelled over the linear space of vertical vector fields 
on the affine jet bundle J l Q — ► i.e., which live in VqJ x Q. Every holonomic connection £ 
defines the corresponding covariant differential on the jet manifold J X Q: 

Dt : J 2 Q-^V Q J l Q C VJ l Q, 

q'oD^O, qioD^qi-C, (7) 
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which takes its values into the vertical tangent bundle VqJ 1 Q of the affine jet bundle J l Q — > 
Q. Any integral section c : R D () — > J l Q for a holonomic connection £ is holonomic, i.e., 
c = c where c is a curve in Q. 

A second order dynamic equation (or simply a dynamic equation) on a configuration 
bundle Q — > R is defined as the kernel, 

qit = e(t,q j ,d), (8) 

of the covariant differential (|7|) for some holonomic connection £ on the jet bundle J X Q — > R. 
Therefore, holonomic connections are also called dynamic equations. By a solution of the 
dynamic equation (§) is meant a curve c in Q whose second order jet prolongation c lives 
in (H). Any integral section c for the holonomic connection £ is the jet prolongation c of a 
solution c of the dynamic equation (§), and wee versa. 

III. NONHOLONOMIC CONSTRAINTS 

Let S be an n-dimensional codistribution on the velocity phase space J x <5- Its annihilator 
Ann (S) is treated as a nonholonomic constraint. Let the codistribution S be locally spanned 
by the 1-forms, 

s a = s^dt + s^dq 1 + %dqi 

on the jet manifold J l Q. Then a dynamic equation £ on the configuration bundle Q — > R is 
said to be compatible with the nonholonomic constraint S if 

s a (O=i\s a = s a + sX + sK l = 0. 

This equation is algebraically solvable for n components of £ iff the n x m matrix s"(g A , ql) 
has everywhere maximal rank n < m. Therefore, we restrict our consideration to the non- 
holonomic constraints, called admissible, such that dimS = dimv(S). 

If a nonholonomic constraint is admissible, there exists a local m x n matrix s l a {q x ,ql) 
such that 

44 = t 

Then the local decomposition ([I]) of a dynamic equation £ can be written in the form 

?=e+ (9) 

The global decomposition ([]]) exists by virtue of the following lemma. 
Lemma 1: The intersection 

W = J 2 Q D Ann (S) 
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is an affine bundle over J X Q, modelled over the vector bundle 

W = VqJ x Q n Ann (S). 

Proof: W consists of the vertical vectors v l d\ G Vq J l Q which fulfill the conditions 

Since the nonholonomic constraint S is admissible, every fiber of W is of dimension m — n, 
i.e., W is a vector bundle, while W is an affine bundle. 

The affine structure of W — > J l Q implies that it has a global section £. 
To construct the global decomposition (TI]), one usually perform a splitting of the vertical 
tangent bundle, 

V Q J l Q = W © V, (10) 
and obtain the corresponding splitting of the second order jet manifold, 

J 2 Q = W © V. (11) 

Here V — > J X Q should be interpreted as the bundle of possible constraint reaction accelera- 
tions. 

If an admissible nonholonomic constraint S is of dimension n = m, a dynamic equation 
£ is decomposed in a unique fashion. If n < m, the decomposition (ffj) is not unique. 
Different variants of this decomposition lead to different constraint reaction forces which, 
from the physical viewpoint, characterize different types of nonholonomic constraints. In next 
Section, we will construct the decomposition of dynamic equations of Newtonian systems 
which corresponds to ideal nonholonomic constraints. 

Now, let us consider some important examples of nonholonomic constraints. 

Let N be a closed imbedded submanifold of the velocity phase space J X Q, defined locally 
by the equations 

f a (q\q$)=0, a = l,...,n<m. 

One can treat iV as a nonholonomic constraint given by the codistribution S = Ann (TN) 
on J X Q | at. This codistribution is locally spanned by the 1- forms 

s a = df a = d t f a dt + diPdqi + d]f a dqi. 

The nonholonomic constraint iV is admissible iff the matrix (<9*/ a ) is of maximal rank n. It 
follows that N is a fibred submanifold of the affine jet bundle J X Q — > Q. 



5 



A nonholonomic constraint N is said to be linear if it is an affine sub-bundle of the affine 
jet bundle J l Q Q. Locally, a linear constraint TV is given by the equations 

f a = / V) + fttfH = °> (12) 

where the matrix /" is of maximal rank. A linear constraint is always admissible. Since iV 
is an affine sub-bundle of J X Q — > Q, it has a global section T which is a connection on the 
configuration bundle Q — > R, called the constraint reference frame. Then, the constraint 
equations ([12]) take the form 

f?(q X )(ql-r)=0. (13) 

We can say that the linear constraint is immovable with respect to the constraint reference 
frame T. Then, one can think of q l T = ql — T l , satisfying the equation flTBP , as virtual velocities 
relative to the linear constraint N. 

Let now a configuration space Q admit a composite fibration Q —>■ £ — > R, where 

ttqs : Q — > S 

is a fiber bundle, and let (t,a r ,q a ) be coordinates on Q, compatible with this fibration. 
Given a connection, 

B = dt ® (dt + B a d a ) + dcr r ® (d r + B?d a ), (14) 

on the fiber bundle Q — > S, we have the corresponding horizontal splitting of the tangent 
bundle TQ. Restricted to the jet manifold J l Q C TQ, this splitting reads 

J X Q = ^(Tr^J 1 ^)©^^, 
Q 

a 4 + a[9 r + % ffl 9 a = + 5 a a a ) + a r t {d r + B«9 a )] + [q a t - B a - a r t B?]d a , 

where tiq^J 1 ^ is the pull-back of the affine jet bundle J 1 !! — > S onto It is readily 
observed that 

iV = B(7T* Q x J l Z) 

is an affine sub-bundle of the affine jet bundle J l Q — * Q, defined locally by the equations 

q1-alB a r {q X )-B a {q x ) = Q. 

This sub-bundle yields a linear nonholonomic constraint. 16,17 The corresponding codistribu- 
tion S = Ann (TN) is locally spanned by the 1-forms, 

s a = -(d t B a + a r t d t B a r )dt - {d s B a + a 7 t d s B a r )da s - (15) 
{d b B a + a r t d b B a r )dq h + dq a t - B a r da r t . 
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With the connection (H), we also have the splitting of the vertical tangent bundle VQ of 
Q — > R and the corresponding splitting of the vertical tangent bundle VqJ 1 Q which reads 



VqJ'Q = W © V, 

^ + g t a a* = of (flf + 5 P fl #) + (# - B?a r t )dl. (16) 

It is readily observed that W |jy consists of vertical vectors which are the annihilators of the 
codistribution (|1|). The splitting fllBD yields the corresponding splitting ([□]) of the second 
order jet manifold J 2 Q. Then we obtain the decomposition ([]]) of every dynamic equation £ 
on J X Q as 

r = r, r = e a - AO- 

IV. NEWTONIAN SYSTEMS WITH NONHOLONOMIC CONSTRAINTS 

Let Q — > R be a fiber bundle together with (i) a non-degenerate fiber metric, 

m: J l Q^ V*Q <g> F*Q, m = -mydg* V dg J ', 

in the fiber bundle Vq J x Q — > J X Q which satisfies the symmetry condition, 

dlrriij = djrriik, (17) 

and (ii) a dynamic equation £ @ on the jet bundle J 1 ^ — ► R, related to the fiber metric m 
by the compatibility condition, 

2£J dm, + m lk d]e + m ife d*£ fc = 0. (18) 

The triple (Q,m,£) is called a Newtonian system. 15 A Newtonian system is said to be stan- 
dard if m is the pull-back of a fiber metric in the vertical tangent bundle VQ in accordance 
with the isomorphism (f|). In this case, m is independent of the velocity coordinates q\. 

The notion of a Newtonian system generalizes the second Newton law of particle me- 
chanics. Indeed, the dynamic equation for a Newtonian system is equivalent to the equation 

m lk (qt t -e)=0. (19) 

Therefore, m is called a mass metric. 

There are two main reasons in order to consider Newtonian systems. From the physical 
viewpoint, with a mass metric, we can introduce the notion of an external force, defined as 
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a section of the vertical cotangent bundle VqJ 1 Q — > ^Q. Let (Q,m,£) be a Newtonian 
system and F an external force. Then 

e F = e + (m-^Fk, 

is a dynamic equation, but the triple (Q, m, £p) is a Newtonian system only if F possesses 
the property 

d\F 3 + d]Fi = 0. (20) 

From the mathematical viewpoint, the equation (|T9| ) is the kernel of an Euler-Lagrange-type 
operator. By an appropriate choice of a mass metric, one may hope to bring it into Lagrange 
equations. This is the well-known inverse problem in time-dependent mechanics. 

Here, we consider Newtonian systems because they provide the vertical tangent bundle 
Vq J l Q with a nondegenerate fiber metric m. Let us assume that m is a Riemannian metric. 
With this metric, we immediately obtain the splitting (|10| ), where V is the orthocomplement 
of W. Then the corresponding decomposition (||) takes the form, 9 

e = e+m ab m^sp b (o, (2i) 

where m a b is the inverse matrix of 

m ab = sfs b ,m ij . 



It is readily observed that the decomposition (piD satisfies the generalized D'Alembert 
principle. The constraint reaction acceleration, 

-r* = -m afe m^s fe (0, (22) 

is orthogonal to every element of Vq J x Q H Ann (S) with respect to the mass metric m. 
Since elements of VqJ 1 QC\ Ann (S) can be treated as the virtual accelerations relative to the 
nonholonomic constraint S, the constraint reaction acceleration ( p2|) characterizes S as an 
ideal constraint. 

The Gauss principle is also fulfilled as follows. Given a dynamic equation £ and the 
above-mentioned fiber metric m, let us define a positive function G(w) on J 2 Q as 

G(w) = m (f (7Ti(w)) - w, f (7r?(w)) - io) , 

G(q\ qt, q\ t ) = m^, % fc )(e(g A , q k t ) - > q k t ) - it)- 

We say that ||w|| = G(w) l l 2 is a norm of w G J 2 Q. 

Proposition 2: Among all dynamic equations compatible with a nonholonomic constraint, 
the dynamic equation £ defined by the decomposition (pip is that of least norm. 
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Proof: Let £ be another dynamic equation which takes its values into W. Then £ — ( C W 
and 

m(f-c,e-e) = o. 

Hence, we obtain 

iicii = m(e-f+f-c,e-f+f-c) = lien + m(f - c, f - o- 

In next Section, we will show that, in the case of nondegenerate Lagrangian systems and 
linear nonholonomic constraints, the decomposition ( pT|) satisfies the traditional D'Alembert 
principle. 

IV. LAGRANGIAN SYSTEMS WITH NONHOLONOMIC CONSTRAINTS 

Nondegenerate Lagrangian systems are particular Newtonian systems. 
A Lagrangian is defined horizontal density, 

L = Colt, C:J l Q^ R, (23) 

on the velocity phase space J l Q. Here, we apply in a straightforward manner the first 
variational formula. 13 ' 15 

Let us consider a projectable vector field 

u = u l d t + u l di, u* = 0, 1, 

on the configuration bundle Q — > R and calculate the Lie derivative of the Lagrangian (|23|) 
along the jet prolongation, 

u = u l dt + rfdi + dttfdj, 

of u, where dt = dt + q\di + • • • is the operator of formal derivative. We obtain 

LuL = (u\ dC)dt = (u*dt + u% + dtu^DCdt. (24) 

The first variational formula provides the following canonical decomposition of the Lie 
derivative fl24|) in accordance with the variational problem: 

u\ dC = {u { - vtqfiSi + d t (u\ H L ) (25) 

where 

H L = v(dL) + L = TTidq* - (tt^J - C)dt (26) 
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is the Poincare-Cartan form and 



£ L : J 2 Q -> V*Q, 

S L = £.p = (d t - d t dt)LB\ (27) 
is the Euler-Lagrange operator for L. We will use the notation 

iii = d^C, iiji = d)d\L. 
A Lagrangian L is called nondegenerate if det rtji 7^ everywhere on the velocity phase space 

•/'a 

The kernel KerSi C J 2 Q of the Euler-Lagrange operator Q27D defines the system of 
second order differential equations, 

ft - d t df)C = 0, (28) 

on Q, called the Lagrange equations. Their solutions are (local) section c of the fiber bundle 
Q — » R whose second order jet prolongations c live in (|28|). 

A holonomic connection on the jet bundle J l Q — > R is said to be a Lagrangian connection 
£l for the Lagrangian L if it takes its values in the kernel (|28|) of the Euler-Lagrange operator 
Sl- Every Lagrangian connection £x defines a dynamic equation on the configuration space 
Q whose solutions are also solutions of the Lagrange equations ( p8|) . If L is non-degenerate, 
the Lagrange equation ( p8|) can be algebraically solved for the second order derivatives, and 
they are equivalent to the dynamic equation, 

= & a = (T" 1 )% + 3«. (29) 

called the Lagrange dynamic equation. 

Every Lagrangian L on the jet manifold J X Q yields the Legendre map, 

L:J l Q^ V*Q, Pi oL = 7Tj, (30) 

where {t,q l ,pi) are holonomic coordinates on the vertical cotangent bundle V*Q. As is well 
known, the Legendre map ( |30D is a local diffeomorphism iff L is nondegenerate. A Lagrangian 
L is called hyperregular if the Legendre map L is a diffeomorphism. 
The vertical tangent map VL to the Legendre map L reads 

vl : Vn j x <5 -> w*Q s v*g x v*g. 

Q 

It yields the linear fibred morphism VqJ 1 ^ — > VqJ l Q and the corresponding mapping, 

J'Q -> ^J X Q ® V3J X Q, my = Try. (31) 
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If a Lagrangian L is nondegenerate, then ( j31"D is a mass metric which satisfies the symmetry 
condition (|17D and the compatibility condition ( |l"8D for the Lagrange dynamic equation ( |29|) . 



Thus, every nondegenerate Lagrangian L defines a Newtonian system. Moreover, a non- 
degenerate Lagrangian system plus an external force which fulfills the condition (|20|) is also 
a Newtonian system. Conversely, every standard Newtonian system can be seen as a La- 
grangian system with the Lagrangian, 

L=^ miJ (ql-r)(qi-P)dt, (32) 

where T is a reference frame, plus an external force. 

Given a nondegenerate Lagrangian L with a Riemannian mass metric = 7Tjj, let now 
S be an admissible nonholonomic constraint on the velocity phase space J X Q. Since this is 
a particular Newtonian system, we obtain the dynamic equation 

Qlt = ?i - m ab m^ s^sUl + « + 4), (33) 
(i^m^i-dt^-dkTT^ + djjr,), 

which is compatible with the constraint S, treated as an ideal nonholonomic constraint. This 
is the Lagrange dynamic equation in the presence of the additional constraint reaction force 

F{ = -m ab sts\£ L ). (34) 

Let us consider the energy conservation law in the presence of this force. 

The energy conservation law in Lagrangian time-dependent mechanics is deduced from 
the first variational formula (^) when the vector field u = T is a reference frame. On the 
shell £i = (|28|), this formula leads to the weak identity, 



LpL w —dt^iTiq^ — £), (35) 

where q^ = q\ — T l is a relative velocity and 

T r = TTiQf, - C (36) 

is the energy function with respect to the reference frame r. 14 ' 15 ' 18 In the presence of an 
external force F, i.e., on the shell Si = — F i: the weak identity (R5j) is modified as 



LpL - q^Fi = -d t T r . 

It is readily observed that, if a nonholonomic constraint is linear and T is a constraint refer- 
ence frame, the constraint reaction force (|34"1) does not contribute to the energy conservation 
law. It follows that, in this case, the standard D'Alembert principle holds, while the equation 
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(|33| ) describes a motion in the presence of an ideal nonholonomic constraint in the spirit of 
this principle. 

The constrained equation of motion (|33| ) is neither Lagrange equations nor a dynamic 
equation of a Newtonian system. In Section VI, we aim to show that it can be seen as a 
part of Hamilton equations in the framework of the Hamiltonian formalism extended to the 
configuration space VQ. 

V. VERTICAL EXTENSION OF HAMILTONIAN FORMALISM 

This Section provides a brief exposition of Hamiltonian formalism of time- dependent 
mechanics on a configuration bundle Q — > R and its extension to the vertical configuration 
space VQ. We consider this extension because any first order dynamic equation on the 
momentum phase space V*Q can be seen as a Hamilton equation in the framework of the 
extended Hamiltonian formalism. This extension is also of interest in the path-integral 
formulation of mechanics. 19 ' 20 

Given a mechanical system on a configuration bundle Q — > R, its momentum phase space 
is the vertical cotangent bundle V*Q of Q — > R, equipped with the holonomic coordinates 
(t, q l ,Pi = (ji). 13-15 The momentum phase space V*Q is endowed with the canonical exterior 
3-form, 

Q = dpi A dq l A dt. 

Let us consider the cotangent bundle T*Q of Q with the holonomic coordinates (t, q l ,p,pi). 
It admits the canonical Liouville form 

E = pdt + pidq\ (37) 

An exterior 1-form H on the momentum phase space V*Q is called a Hamiltonian form if it 
is the pull-back 

H — h*E — pidq 1 - Hdt (38) 
of the Liouville form H (|3T| ) by a section h of the fiber bundle 

C : T*Q -> V*Q. (39) 

Remark: With respect to a trivialization Q = R x M, the Hamiltonian form ( |3"%| ) is the well- 
known integral invariant of Poincare-Cartan, where Ti. is a Hamiltonian. The peculiarity of 
Hamiltonian time-dependent mechanics issues from the fact that Hamiltonians are not scalar 
functions under time-dependent transformations, but make up an affine space modelled over 
the linear space of functions on V*Q. 
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For instance, every connection T on a configuration bundle Q — > R is an affine section, 

poT = -piT\ 

of the fiber bundle fl39|), and defines the Hamiltonian form 

Hr = Pidq 1 - piVdt. 

It follows that any Hamiltonian form on the momentum phase space V*Q admits the split- 
ting, 

H = H r - H v dt = pidq 1 - (^P + H r )dt, 

where T is a connection on Q — ► R and 7Yr is a real function on V*Q, called the Hamiltonian 
function. The following assertions are basic facts in the Hamiltonian formulation of time- 
dependent mechanics. 14,15 

Proposition 3: Every Hamiltonian form H on the momentum phase space V*Q defines the 
associated Hamiltonian map, 

H : V*Q -> J X Q, (f t oH = Shi. 



Proposition 4-' Given a Hamiltonian form H on the momentum phase space V*Q, there 
exists a unique connection 

j H = dt + d l Hdi - diHd 1 (40) 
on V*Q — > R, called a Hamiltonian connection, such that 

7h J n = dH. 

The kernel of the covariant differential of the Hamiltonian connection ([IDp defines the 
Hamilton equations, 

q\ = &"H, (41a) 
Pu = -W, (41b) 

for the Hamiltonian form H. Their solutions are integral curves for the Hamiltonian con- 
nection j H pop. 
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Now let us consider the vertical tangent bundle VQ of the fiber bundle Q — > R, coordi- 
nated by (t,q\q l ). It can be seen as a new configuration space, called the vertical configu- 
ration space. The corresponding vertical momentum phase space is the vertical cotangent 
bundle V*VQ of VQ — > R. The vertical momentum phase space V*VQ is canonically 
isomorphic to the vertical tangent bundle VV*Q of the ordinary momentum phase space 
V*Q — > R, coordinated by (t,q t ,Pi,q l ,pi). It is easily seen from the transformation laws 
that (q l ,pi) and (q\Pi) are canonically conjugate pairs. 

The vertical momentum phase space VV*Q is endowed with the canonical 3-form, 

Oy = [dpi A dq l + dpi A dq l ] A dt. 

For the sake of brevity, one can write f2y = <9yf2, where <9y = <f <9j + pj<9* is the vertical 
derivative. 

The notions of a Hamiltonian connection, a Hamiltonian form, etc., on the vertical mo- 
mentum phase space VV*Q = V*VQ are introduced similarly to those on the ordinary 
momentum phase space V*Q. In particular, a Hamiltonian form on VV*Q reads 

H v = pidq 1 + pidq 1 — Jiydt. 

Since Hamiltonian forms are determined modulo exact forms and the function pitf is globally 
defined on VV*Q, we will write 

H v = Pidq 1 - q l dpi - Hydt. (42) 
The corresponding Hamilton equations read 





= q\ = d'Hv, 


(43a) 


ii 


= Pu = -diH v , 


(43b) 


¥ 


= ql = d'Hv, 


(43c) 


ii 


= Pa = —diHv, 


(43d) 



where 

7 = ^ + 7^ + 7^ + 7^ + 7^ 

is a Hamiltonian connection on the vertical momentum phase space VV*Q — > R. 

There is the following relation between Hamiltonian formalisms on V*Q and VV*Q. 13 ' 15 
Let VT*Q be the vertical tangent bundle of the cotangent bundle T*Q — >■ R is equipped 
with holonomic coordinates (t,q l ,pi,p,q\pi,p) and endowed with the canonical form, 

E v = pdt + pidq 1 - q l dpi. 
14 



Proposition 5: Let 7# be a Hamiltonian connection on the ordinary momentum phase space 
V*Q — > R for a Hamiltonian form, 



Then the connection 



H = h*E = Pl dq l - Halt. 
Vj H ■ VV*Q -> VJ X V*Q = J 1 VV*Q, 

Vm = d t + fdi + 1% + dvYdi + d vll d\ 



(44) 



(45) 



on the vertical momentum phase space VV*Q — > R is a Hamiltonian connection for the 
Hamiltonian form, 



H v = {Vh)*E v = p t dq l - tfdpi - d v Hdt, 
d v H=(q i d i +p i d i )H, 

where Vh : VV*Q — > VT*Q is the vertical tangent map to h. 
The corresponding Hamilton equations read 



(46) 
(47) 



Y = 


d l H v = 


d l H, 


ii = 


—diHy 


= -dtH, 


r = 


d'Hy = 


d v d l H, 


ii = 


-diHy 


= -dydiU 



(48a) 
(48b) 
(48c) 
(48d) 

It is easily seen that the equations (48a) - ( 48b ) are exactly the Hamilton equations (|41a|) - 
(|41b|) for the Hamiltonian form H. 

Remark: In order to clarify the physical meaning of the Hamilton equations (48c) - ( (48d| ), 
let r(t) be a solution of the Hamilton equations (|48a|) - ( |48b|) . Let f(t) be a Jacobi field, i.e., 
r(t) + er(t) is also a solution of the same Hamilton equations modulo terms of order two in 
e. Then it is readily observed that the Jacobi field r(t) fulfills the Hamilton equations (|48c|) 

-H3). 

The following assertion plays a prominent role in the sequel. 13,15 

Proposition 6: Any connection 7 on the momentum phase space V*Q — ► R gives rise to the 
Hamiltonian connection, 



Y = l\ li = 7i» f = Pjd'Y - tfd^j, 7^ = -pjda 3 + cpdaj, 
for the Hamiltonian form, 

H v = piidq 1 - fdt) - q\dpi - jidt) = p\dq l - tfdpi - ( Pi f - tf^dt, 



(49) 
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on the vertical momentum phase space VV*Q. 

In particular, if 7 is a Hamiltonian connection on the fiber bundle V*Q — > R, then 
is exactly the connection V7 fl45|) . 

It follows that every first order dynamic equation on the momentum phase space V*Q 
can be seen as the Hamilton equations (|43a|) - (|43b ) for a suitable Hamiltonian form on the 
vertical momentum phase space. 

VI. HAMILTONIAN SYSTEMS WITH NONHOLONOMIC CONSTRAINTS 

Let L be a hyperregular Lagrangian with a Riemannian mass metric m. In this case, 
Hamiltonian and Lagrangian formalisms of time-dependent mechanics are equivalent. There 
exists a unique associated Hamiltonian form H on V*Q such that 

H = L-\ p % = n t (q\m(q\p k )), q\ = d i H{q X , vr i (g A , g t fc )), (50a) 
C o H = 7h J H = p l d l H - H. (50b) 

As an immediate consequence of (|50a|) , we have J 1 if = ( J 1 L) _1 , where the jet prolongations 
of the Hamiltonian and Legendre maps read 

J l H : J l V*Q - J 1 J'Q, (q\ q\, gfo, <&) o J 1 ]? = (g\ TO, g*, d t d i H), 
J 1 ! : JV X Q -> JV*Q, (q X ,Pi,qi,Pti) ° J X L = (g\ tt 4 , gj t) , d^). 

Then, using ( |50a|) - ( |50b| ), we obtain 

7h = J 1 ! o£ L o H. 

Let introduce the notation M iJ = d l cPTi. There are the relations 

rik 1 _ "rj\ ri / 71 /fife 



M l * (m fci o if) = 5} , m fci (M ife o L) — 5], m l3 = n tj . 
It follows that M is a fiber metric in the vertical tangent bundle VqV*Q of the fiber bundle 

Given a codistribution S on J 1 ^, let us consider the pull-back codistribution H*S on 
V*Q, spanned locally by 1-forms 

pa = fl* s a = ( s a + gfify&jQM + ( s « + ^d i d j H)dq i + S^M^dpj = 

ffidt + ffidq 1 + p ai d Pi . 

This codistribution defines a nonholonomic constraint on the momentum phase space V*Q. 
Given a Hamiltonian connection 7# (^), let us find its splitting 

7^ = 7 + ^ (51) 
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where 7 is a connection on V*Q — > R which satisfies the condition 

7CAnn(#*S). (52) 

The connection 7 ( |52"D obviously defines a first order dynamic equation on the momentum 
phase space V*Q which is compatible with the nonholonomic constraint H*S. The decom- 
position (|5"ID is not unique. Let us construct it as follows. 

Given a Hamiltonian connection 7^, we consider the codistribution Sh on V*Q, spanned 
locally by the 1-forms dq % — Y H dt. Its annihilator Ann (Sh) is an affine sub-bundle of the 
affine jet bundle J 1 V*Q — > V*Q, modelled over the vertical tangent bundle VqV*Q. The 
Hamiltonian connection 7# is a section of this sub-bundle. Let us take the intersection 

W = Ann (S H ) H Ann (H*S). 

Lemma 7: W is an affine bundle over V*Q, modelled over the vector bundle, 

W = V Q V*Qf] Ann (H*S). 

Proof: The intersection W consists of elements v = Vid 1 of VqV*Q which fulfill the conditions 

v l (3 ai = 0. 

Since the nonholonomic constraint S is admissible and the matrix M !J is nondegenerate, 
every fiber of W is of dimension m — n, i.e., W is a vector bundle, while W is an affine 
bundle. 

Then, using the fiber metric M in VqV*Q, we obtain the splitting 

V Q V*Q = W®V, 
where V is the orthocomplement of W, and the associated splitting 

Ann (S H ) = W © V. 
The corresponding decomposition ( [51]) reads 

7 = in ~ M ab M l3 P ai (J b ( lH )&, (53) 
where M a b is the inverse matrix of 

M ab = p ai p bj M i:j . 
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The splitting fl53|) is the Hamiltonian counterpart of the splitting (21). We have the 
relations 

~ab = ^ab Q /3 a ( 7H ) = S a (£ L ) O H, 

and as a consequence 

7 = J x LoloH. 

Remark: The above procedure can be extended in a straightforward manner to any standard 



Newtonian system, seen as a Lagrangian system with the Lagrangian (32) and an external 
force. Following this procedure, one may also study a nonholonomic Hamiltonian system, 
without appealing to its Lagrangian counterpart. 

The connection (|53|) defines the system of first order dynamic equations, 

qt = #H, pu = -d-H - M ab M ij p ai b ( r fH), (54) 

on the momentum phase space V*Q, which are not Hamilton equations. Nevertheless, in 
accordance with Proposition |], one can restate the constrained equations of motion fl54|) as 
the Hamilton equations ( |48a| ) - ( [48 b| ) for the Hamiltonian form, 

H v = Pidq 1 - tfdpi - d v Hdt - q i M ab M ij p ai (3 b (j H )dt, 

on the vertical momentum phase space VV*Q, where the last term can be written in brief 
as (-dyJtfJO). 

The Hamiltonian form of the constrained equations of motion may be important in 
connection with the following speculations. 

Given a Hamiltonian form Hy (^) on the vertical momentum phase space VV*Q, let us 
consider the Lagrangian 

Lh = Piq\ ~ tfPu ~ (55) 

on the first order jet manifold J 1 VV*Q of the fiber bundle VV*Q — > R. It is readily observed 
that the corresponding Lagrange equations are exactly the Hamilton equations ( |43a| ) - ( [43 d|) 
for the Hamiltonian form Hy. In particular, let H be a Hamiltonian form on an ordinary 
momentum phase space V*Q and TCy = dyH. In this case, the Lagrangian (|55D reads 

Lh = Pi(qi - &H) ~ q\pu + d{H). 

It is easily seen that this Lagrangian vanishes on solutions of the Hamilton equations for the 
Hamiltonian form H. By this reason, it is applied to the functional integral formulation of 
mechanics. 19 ' 20 
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